Novel Algebraic Boson Liquid phase with soft Graviton excitations 
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A bosonic model on a 3 dimensional fee lattice with emergent low energy excitations, with the 
same polarization and gauge symmetries as gravitons is constructed. The novel phase obtained 
is a stable gapless boson liquid phase, with algebraic boson density correlations. The stability of 
this phase is protected against the instanton effect and superfluidity by self-duality and large gauge 
symmetries. The gapless collective excitation of this phase closely resembles gravitons, although 
they have a soft a; ~ A:^ dispersion relation. The dynamics of this novel phase is described by new 
set of Maxwell equations. This phase also possesses an intricate topological order, requiring 18 
winding numbers to specify each topological sector. 
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The algebraic spin liquid is of crucial importance in 
understanding the various types of behavior of strongly 
correlated electron systems. This is not only because it 
belongs to the family of the spin liquid, which has played 
a central role in the attempt to understand superconduc- 
tivity in the cuprates 1] 2], but also because of its spe- 
cial gapless excitation, which usually resembles one of the 
most fundamental particles in the universe, the photon 
Unlike other types of spin liquids (for instance, 
the Z2 spin liquid of the Quantum Dimer model on a 
triangular lattice 0, and a spin- 1/2 model on Kagome 
lattice j7i]), the algebraic phase behaves like a critical 
point, as the correlation functions of physical quantities 
(for instance, the spin density correlation function) fall 
off algebraically. Therefore when a relevant perturbation 
is turned on, it is supposed to drive the system into an or- 
dered phase. It has been shown that the nature of many 
kinds of normal order can be understood in their mother 
phase, the 2 -fid algebraic spin liquid 8]. 

Since the algebraic spin liquid phase is critical, it is 
usually vulnerable to all kinds of relevant perturbations. 
In 2-|-ld, monopole proliferation gaps the system and 
causes a certain crystalline order to form which breaks 
translation symmetry. Models with QED gauge sym- 
metry can only give rise to gapped gauge bosons 0. 
Monopoles in 2-1- Id can only be irrelevant due to the 
Berry phase at the transition point between different 
types of crystal order in the quantum dimer model 
ISl I 111- or between crystal order and superfluid or- 
der [13, ■ Protecting the algebraic spin liquid phase 
is equivalent to protecting gapless bosonic excitations 
without symmetry breaking. As is well-known, contin- 
uous symmetry breaking is one way to obtain a gapless 
bosonic excitation. The gapless excitation is the Gold- 
stone mode [1^ . However, the algebraic spin liquid which 
is a phase without any symmetry breaking, requires a 
different mechanism to guarantee the stability of its gap- 
lessness and criticality. So far, the only stable algebraic 
spin liquid phase that has been proven to exist is the 
3-f-ld photon spin liquid. This spin liquid can be re- 



alized in either the quantum dimer model on the cubic 
lattice 01, or a spin-1/2 model on the pyrochlore lattice 
4]. In a 2-fld system, monopole proliferation is a con- 
cern, although the existence of a 2-fld stable algebraic 
liquid phase has been argued, by coupling the compact 
gauge fields with a large number of flavors of the matter 
field [3. 

In this paper, a new type of algebraic liquid phase 
is shown to exist. Since a spin system can usually be 
mapped to a bosonic system through the identification. 
Hi — n = sf and fej ^ , the bosonic language is used in 
this paper. A lattice model is constructed, which shows a 
stable algebraic boson liquid phase with gapless collective 
excitations. The gapless collective excitations have the 
same polarization and gauge symmetry as the graviton, 
although the dispersion relation is softened to w ~ k^. 
The stability of this algebraic phase is due to the self- 
duality and large gauge symmetry of the low energy ef- 
fective Hamiltonian. This work focuses on the basic prop- 
erties of this new phase. We leave the discussion of the 2d 
case and the transition to other phases to another piece 
of work ^3 . 

Let us first briefly review some basic properties of the 
3d photon liquid phase as a warmup. Although many 
different models have been proposed in the last few years 
l^^l^, the phases obtained have very similar proper- 
ties. Therefore, let us take a 3d Quantum Dimer Model 
on the cubic lattice as an example. On the cubic lattice, 
every site is shared by six links, and only one of those 
links is occupied, by exactly one dimer. On every square 
face, if two parallel links are occupied by dimers, they 
can be resonated to dimers perpendicular to the origi- 
nal ones. Besides this resonating term, another diagonal 
weight term for each flippable plaquette is also included 
in the Hamiltonian 

H = J2 I + '^•^•) + ^(IIDdll + I =>(= I) (1) 

The dimer model can be mapped onto a rotor model. A 
rotor number can be defined on each link to describe the 
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presence or absence of dimers: n = I when the hnk is oc- 
cupied by a dimer, and n = when the hnk is empty. The 
Hilbert space of this quantum system is a constrained 
one, with the constraint = 1 around each site. 

Let us define the quantity Ei^a as Ei^a = (— l)'^''''""i,a- 
Now, the constraint of this system can be rewritten as the 
Gaussian law for electric fields, = ±1. Notice that 

because the quantity Ei^a is defined on links, a natural 
vector notation can be applied: Ei = {Ei^x,Ei^y,Ei^z)- 

The background charge distribution ±1 on the cubic 
lattice plays a very important role in the solid phase, i.e. 
the confined phase. The form of the crystalline phase can 
be determined from the Berry phase induced by the back- 
ground charge distribution |l3l||3|. However, as we are 
focusing on the algebraic liquid phase, the background 
charge is not very important. Let us instead impose the 
constraint d^E^ = 0. Because of this local constraint, 
the low energy physics will be invariant under the gauge 
transformation A — > A + V f {A is the conjugate variable 
of E), which is exactly the gauge symmetry of the U{1) 
gauge theory. Now the effective Hamiltonian of this the- 
ory should be invariant under this gauge transformation. 
The Hamiltonian can be effectively written as 

H = ^-icos(V X 1) 1/(2k)£;2 (2) 

This Hamiltonian is the 3-|-ld compact QED, which 
should have a deconfined photon phase llg)- The 
reason for the existence of the photon phase is actually 
the self-duality and gauge symmetry, as discussed below. 

First, because the constraint d^E^ = is strictly 
imposed on this system, the matter field, i.e. the de- 
fect which violates this constraint, is absent. As is well 
known, the U(l) gauge symmetry cannot be broken spon- 
taneously without a matter field Therefore the 
dimer superfluid order is ruled out. The other instability 
of the photon phase is towards the gapped solid phase, 
which can be analyzed in the dual theory. We can in- 
troduce the dual vector h and dual momentum vector tt 
[h and tt are both defined on the faces of the cubic lat- 
tice) as E = V X h and V x A = 7?. One can check the 
commutation relation and see that h and tt are a pair of 
conjugate variables. The Gaussian law constraint on this 
system is automatically solved by the vector h. Now, the 
field theory for the photon phase is self-dual: 

H = -ti{V X A)^ + ciE^y -E^O, 

= -t2(V X +C27f^V • 7? = 0. (3) 

In principle, a vertex operator cos{2TTNh) is supposed to 
exist in the dual Hamiltonian, due to the fact that E only 
takes on integer values. N is an integer corresponding to 
the Berry phase of the vertex operator. When this vertex 
operator is relevant, it will gap out the photon excitation 
and drive the system into a crystalline phase, according 
to the Berry phase. However, the vertex operator is ir- 
relevant in the photon phase. Notice that, because the 
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FIG. 1: The xy plane of the fee lattice. On eaeh site there 
are three orbital levels, we denote the particle numbers as 
(ni,n2,n3); on each face center there is one orbital level, we 
denote the particle number as n. 

theory is self-dual, the dual theory has the same gauge 
invariance h h + V f as the original theory. However, 
the vertex operator cos{2nNh) breaks gauge invariance 
and thus the correlation function between two vertex op- 
erators in this photon phase is zero, i.e. the vertex op- 
erator is irrelevant in this Gaussian theory. Unless the 
dual charges (the monopole) condense and break the dual 
U{1) gauge symmetry, the photon phase is always stable. 
The gaplessness of this phase is protected by both self- 
duality and the gauge symmetry. 

The photon phase is an algebraic liquid phase, the 
dimer density-density correlation function falls off alge- 
braically. When t = V in the lattice model can 
be solved exactly, and the equal-time density correlation 
falls off as (n(O)n(r)) - l/r^ 

Let us now move on to our graviton model. A 3d fee 
lattice is considered and physical quantities are defined 
on both the sites and face centers of the cubic lattice. Let 
us assume there are 3 orbital levels on each site, and 1 
orbital level on each face center. The particle number on 
the face center is denoted as n, and the particle numbers 
on sites are denoted as (ni, 712, n^) in Fig^ 

The Hamiltonian for this system contains three parts, 
H = Hq + Hi + H2- Hi is the nearest neighbor 
hopping between sites and their nearest face centers, 
and also between adjacent face centers (notice that 
adjacent face centers have the same distance between 
them as a site and its nearest face center). Hi = 

(E<.J> ELi ~tbl^b] - E<u> + h-^-)^ is the 
on site interaction H2 — U{n — 1)^, which fixes the av- 
erage filling of the fee lattice. Hq is the interaction term 
involving links in all 3 directions. For example, for the 
link in («,a;), the interaction term reads 

Hq = ^(f^i+l/2^ + l/2z + '^i+l/2x-l/2z + ^i+l/2£+l/2y 

+"i+i/2£-i/2{i + 2ni,i + 2ni,i+£ - 8)^ (4) 
. The links in y and z directions are treated similarly. If 
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FIG. 2: The distribution of the sign 77 on xy plane. After 
introducing 77, the constraint can be written as in (|SJ. 
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The cosine terms involving the curvature tensor are ring 
exchange terms generated by the nearest neighbor hop- 
ping, resembling the dimer flipping term cos(V x A) in 
(0. Ring exchanges in © are generated at the eighth 
order perturbation of the nearest neighbor hopping, ii, 
^2 ~ jV'^ . All lower order perturbations only generate 
terms which do not comply with the constraint One 
of the ring exchanges is shown in Fig. ^ it corre- 
sponds to the term cos{Rxyxy) ■ 

In order to derive the correct action, one needs to in- 
troduce a Lagrange multiplier AiQ for the constraint jSjl . 
The full action reads 



the bracket in Q is expanded, it becomes the usual two 
body repulsion term. 

When Hq becomes the dominant term in the Hamil- 
tonian, it effectively imposes a constraint on the system. 
The best way to view this constraint is by introducing a 
staggered sign and defining new variables, similar to the 
electric field in the dimer model discussed before. Let us 
define a rank-2 tensor Eij. The off-diagonal terms are 
defined on face centers as Eij — rir{n — 1). n is located 
at one of ij face centers; the diagonal term is defined on 
sites as En = "qri^i — 1) with i = 1,2,3. 77^ = ±1 and 
the distribution of sign r/^ is shown in Fig. |21 

After introducing the sign ry, the constraint effectively 
imposed by can be compactly written as 
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A violation of this constraint can be interpreted as a 
charged defect excitation and can drive the system into 
an ordered boson superfluid state after condensation. We 
will discuss this transition in another paper [l^ . In 
the current work we only focus on the case when the 
constraint is strictly imposed. The constraint ^ re- 
quires the low energy Hamiltonian to be invariant un- 



der the gauge transformation 



Aij + difj + djfi 



. This is precisely the gauge symmetry of the graviton 
in 3d space. Thus, the low energy physics can only in- 
volve the gauge invariant quantity Ra^f3i, — ^/'2,{Aa^^^i3 + 



L = ^{drA^j - diAjo - drAiof + ^ ii cos(i?yij) 

+ ^ i2Cos(i?yife)(7) 

Now, the gauge symmetry can be enlarged to quantities 



defined in space-time: A^ 



A 



d/j-ff + dufti and 



^00 = 0, /o = 0. In this system, the boson superfluid or- 
der is again ruled out by the gauge symmetry. Without 
crystalline order (proven later), the system is in a liquid 
phase with excitations which have the same gauge sym- 
metry (and hence the same polarization) as the graviton. 
Unlike the quantum dimer model, the curvature tensor is 
a second spatial derivative of Aij and the graviton mode 
has a soft dispersion a; ~ fc^. 

Whether the graviton excitations survive, or crystal 
order develops can be determined in the dual theory. 
(5) If we define the symmetric tensor £ij as En = V^Ei,, 



£ij = l/\/2Eij,i ^ j, the constraint ISJ can be solved 
by defining the dual tensor hij as £ij — eiab^jcddadchbd- 
This is a double curl of the symmetric tensor hij. hij 
also lives on the sites and faces of this fee lattice. 

If checked carefully, one can notice that, the curvature 
tensor can also be written in the double curl form 



2Rxyxy ^zab^zcd^a^cAltd^ 2,Rx 



2R. 



ryxy 
yzyz 



^yab^ycd^a ^c-^i 



~ ^xab^xcd^adc-^hdi'^Rxyxz — ^yah^zcd^a^c-^\ 



'^Ryxyz — ^xab^zcd^a^c-^bdj '^-^ 



■xyxz 
zxzy 



bd-> 
^bd-> 



xab^ycd 



Therefore, this model is precisely self-dual, as long as 



^0/3,^1/ ), the linearized curvature ten- we define the dual variables hi 



sor. Here the convention of general relativity has been 
used: A^^^^^p = df^dpA^^. 

In 3d space, the curvature tensor has 6 nonzero com- 
ponents, because the total number of nonzero compo- 
nents is reduced by the symmetry of curvature tensor 



in terms of £ij ~ 



f-iabf-jcddadjibd and its conjugate Hij as follows: 

'^/^Rxyxy — '^zz^ \'^Ryzyz — '^xx^ \'^Rxzxz — '^yy 
lyf^Rxzyz — "^xy^ 2^\/2Rxyxz — '^yz^ 2y/2Rxyzy — "^xz 



(9) 
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Thus, now the T^ij is subject to the same constraint as 



effective low energy Hamiltonian reads 



Hring ^ ^ ^1 ^^^{Rijij 



^ i2Cos{Rijik) 

ijk,i^jj^k,i^k 



Now the dual action reads 

Lduai — y^{9thij — dihjo — djhio)'^ — ^ piR 



ijij 
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•(10) 



RijU is the curvature tensor of hij . hio is a Lagrange mul- 
tiplier, due to the constraint on TTy . The ellipses include 
possible vertex operators. Without the vertex operators, 
this theory is at a Gaussian fixed point and hence in an 
algebraic liquid phase with soft graviton excitations. If 
the vertex operators are relevant, they will destabilize 
the liquid phase and gap the graviton excitation, and 
form crystalline order according to the Berry phase. The 
dual action I|1U|I is also invariant under the gauge trans- 
formation hfj,i, h^i, + d^f^ + d^ffj,, with hoo = and 
fo = 0. Therefore, any kind of vertex operator (for exam- 
ple cos{2NTrhij)) is not a gauge invariant operator. The 
correlation function between two vertex operators at the 
Gaussian fixed point is zero or correlated at very short 
range, so the Gaussian fixed point (also the algebraic spin 
liquid) is stable against weak perturbations of the vertex 
operators. 

In the critical liquid phase, correlation functions be- 
tween local order parameters decay algebraically. In our 
case, the correlation function decays more rapidly than 
in the QED model, because the local order parameter is 
the second order derivative of the dual variable hij . For 
instance, the correlation between two spatially distant 
boson density operator ni—n scales as ~ (— l)^+^+^/r^. 

The dynamics of the liquid phase can be described by 
a new set of Maxwell equations. Define the rank-2 tensor 
= tiab^jcddadcAbd- The dynamical equations that 
describe this liquid phase are 



F) F — (9 K 

K^iab^icddadcBbd 



dtBij + K€iab£jcddadc£bd 



0, 
0, 
0. 



(11) 



Charged excitations and topological defects are absent 
in these equations, so they correspond to the Maxwell's 
equations in vacuum. If the constraint © is softened, 
charge density and charge current have to be incorpo- 
rated in equation Because the field variables are 
rank-2 tensors, the charge density will be a vector field, 
and the charge current will again be a rank-2 tensor. 
Therefore, the form of charge density and current are 
very similar to spin density and spin current. We will 
leave the detailed discussion to the other paper 0| . 

Since the spin liquid does not break any symmetry, 
one cannot classify spin liquids by symmetry. Instead, 
topological order has been introduced to classify differ- 
ent spin liquids jij. The topological sector of the spin 
liquid phase with photon excitations can be described by 
6 integers, 3 of them describe electric flux through the 
system and the other 3 describe the magnetic flux ^ . In 
our case there are more topological sectors. Consider a 
yz plane in the original lattice and sum over E^x of a 
particular configuration on this plane. This gives us an 
integer. This integer is invariant under any ring exchange 



allowed. In addition, this quantity does not depend on 
which yz plane is chosen, this is due to the constraint 
© and Gauss's theorem, / (PxdiCix = J Cix ■ dSi = 0. 
The total integral of Ei^ on any yz plane gives the same 
result. Let us denote this integral as nE,xx,yz- 

There are similar integrals which also specify a topolog- 
ical sector: they are nE,yy,xz, ns^zz^xy, nE,xy,yz, ns.xy.xz, 
nE,yz,xz, nE.yz.xy, nE,xz,xy, riE,xz,yz- There are in total 9 
types of electric fiux winding numbers. The dual field TTy 
is subject to the same type of constraint as ©, because 
of the self-duality. Thus, the flux of TTy provides another 
9 integrals specifying a topological sector. In total there 
are 18 winding numbers. 

To conclude, in this work a new type of stable alge- 
braic boson liquid phase has been realized in a 3d lattice 
model. Low temperature physics should be controlled by 
the gapless graviton excitations, for instance, the gravi- 
tons make the contribution to the specific heat C ^ T'^/^, 
which is larger than the T'^ contribution from phonons. 
If the original model was written in terms of fermionic 
particles, novel non-Fermi liquid behavior is expected. 

The author would like to thank L. Balents, M. P. A. 
Fisher, A. Vishwanath, J. E. Moore and C. Wu for helpful 
discussions. 
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